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After the early suggestion by John Pendry to probe unoccupied bands at surfaces through the 
time reversal of the photoemission process, the inverse-photoemission technique yielded the first 
conclusive experimental evidence for the existence of image-potential bound states at metal surfaces 
and has led over the last two decades to an active area of research in condensed-matter and surface 
physics. Here we describe the current status of the many-body theory of inelastic lifetimes of 
these image-potential states and also the Shockley surface states that exist near the Fermi level in 
the projected bulk band gap of simple and noble metals. New calculations of the self-energy and 
lifetime of surface states on Au surfaces are presented as well, by using the GWT approximation of 
many-body theory. 

PACS numbers: 71.45.Gm, 78.68.+m, 78.70.-g 



I. INTRODUCTION 

In a pioneering paper [![, Echenique and Pendry in- 
vestigated the observability of Rydberg-like electronic 
states trapped at metal surfaces via low-energy electron 
diffraction (LEED) experiments. They discussed the life- 
time broadening of these image-potential-induced surface 
states (image states), and reached the important con- 
clusion that they could, in principle, be resolved for all 
members of the Rydberg series. 

A few years later, Pendry suggested a new experi- 
ment measurement of the bremsstrahlung-radiation 
spectrum from electrons, with energies no more than 
a few tens of electron volts, incident on clean surfaces, 
thereby turning incident electrons into emitted photons. 
This photon-emission experiment is simply the time re- 
versal of the photoemission process and was referred to 
by Pendry as inverse photoemission, or IPE for short. 

Subsequently, Johnson and Smith 3] pointed out 
that image states were potentially observable by angle- 
resolved IPE; using this technique, Dose et al. [4| and 
Straub and Himpsel [|[ reported the first conclusive ex- 
perimental evidence for image-potential bound states at 
the (100) surfaces of copper and gold. Since then, several 
observations of ima ge s tates have been made using this 
technique [1, @, H, foL IXQl, [TTI . Il2| , and also the more recent 
high-resolution techniques of two-photon photoemission 
(2PPE) [H, Q [ll| and time-resolved two-photon pho- 
toemission (TR-2PPE) [U, EE El- In 2PPE, intense 
laser radiation is used to populate an unoccupied state 
with the first photon and to photoionize from the inter- 
mediate state with the second photon. In TR-2PPE, the 
probe pulse which ionizes the intermediate state is de- 
layed with respect to the pump pulse which populates it, 
thus providing a direct measurement of the intermediate- 
state lifetime. 

At metal surfaces, in addition to image states (which 
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FIG. 1: The TL projected bulk band structure (shaded areas) 
of the (111) and (100) surfaces of the noble metals Cu and 
Au. The solid lines represent Shockley (n = 0) and image- 
potential (n = 1) surface-state bands. Well-defined Shock- 
ley states are only present at the Cu(If I) and Au(lll) sur- 
faces, where the projected band gap extends below the Fermi 
level. Well-defined image states are present at the Cu(f00) 
and Au(f 00) surfaces, and also at the (111) surface of Cu. At 
Au(lll) there is only an image-state resonance lying above 
the top of the band gap. The horizontal thin solid lines rep- 
resent Fermi and vacuum levels. 



are originated in the combination of the long-range im- 
age potential in front of solid surfaces with the pres- 
ence of a band gap near the vacuum level) [H, [2(| 
there exist crystal-induced surface states (which would 
occur even for step barrier in the absence of the im- 
age potential) 



2l| often classified as Shockley and Tamm 



states [22], [23j : Shockley states exist near the Fermi level 
in the projected bulk band gap of simple and noble met- 
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als, and Tamm states exist at the M points of the surface 
Brillouin zone for various noble- metal surfaces. The life- 
times of excited holes at the band edge (k|| = 0) of Shock- 
ley states have been investigated with hi gh-r esolution 
angle-resolved photoemission (ARP) [U [H, [H, H3] 
and with the use of the scanning tunneling microscope 
(STM) [H, [li]. STM techniques have also allowed the 
determination of the lifetimes of excited Shockley and 
image electrons over a range of energies above the Fermi 
level [HEU. 

Figure Q] illustrates Shockley and image-potential 
states in the gap of the TL projected band structure of 
the (100) and (111) surfaces of the noble metals Cu and 
Au. If an electron or hole is added to the solid at one of 
these states, inelastic coupling of the excited quasiparti- 
cle with the crystal may occur through electron-electron 
(e-e) and electron-phonon (e-ph) scattering. 

In this paper, we first give a brief description of exist- 
ing calculations of e-ph inelastic linewidths of image and 
Shockley states, and we then focus on the many-body 
theory of e-e inelastic lifetimes of these states. In partic- 
ular, we describe the current status of many-body GW 
and GWT calculations, and we report new GWT calcu- 
lations of the self-energy and lifetime of surface states 
on Au surfaces. We conclude that short-range exchange- 
correlation (xc) contributions to the electron (or hole) 
self-energy are small, as occurs in the case of bulk states. 

Unless otherwise is stated, atomic units are used 
throughout, i. e., e 2 — h — m e = 1. The atomic unit 
of length is the Bohr radius, ao = h 2 /to 2 = 0.529A, the 
atomic unit of energy is the Hartree, 1 Hartree = e 2 /ao = 
27.2 eV, and the atomic unit of velocity is the Bohr ve- 
locity, Vo — ac — 2.19 x 10 8 cms~ 1 , a and c being the 
fine-structure constant and the velocity of light, respec- 
tively. 



II. ELECTRON-PHONON COUPLING 

The decay rate due to the e-ph interaction, which is 
relatively important only in the case of excited Shock- 
ley holes near the Fermi level, has been investigated 
by using the Eliashberg function [32| . In particular, at 
zero temperature (T = 0) and in the high-temperature 
limit (ksT » uj m , ks being Boltzmann's constant and 
Lu rn , the maximum phonon frequency) one finds respec- 
tively, assuming translational invariance in the plane of 
the surface, the following expressions for the e-ph induced 
linewidth (or lifetime broadening) T ep of surface states of 
parallel momentum k and energy E [32l |33| : 

/■max(|E| ,uj m ) 

T° ep (k 7 E) = 2vr / a 2 F^{u)duj (1) 

Jo 

and 

r ep (k,£;) = 2n\(k)k B T, (2) 



TABLE I: Electron-phonon linewidths at T = {T%), in 
meV, and the parameter A that is responsible for the high- 
temperature behaviour of T ep , as reported at the T point in 

Refs. mm- 



Surface n F° p A 

Al(100) 18 023 

Cu(lll) 7.3 0.16 

Ag(lll) 3.7 0.12 

Au(lll) 3.6 0.11 

Cu(100) 1 <lmeV < 0.01 

Ag(100) 1 <lmeV < 0.01 



where a 2 Fk(w) is the Eliashberg function, which repre- 
sents a weighted phonon density of states, and 

Jo u 

For many years, the e-ph contribution r ep (k, E) to the 
inelastic decay of surface states had been calculated using 
a three-dimensional (3D) Debye phonon model with A ob- 
tained from measurements or calculations of bulk prop- 
erties [32j. More refined calculations, which are based 
on an accurate description of the full Eliashberg spectral 
function, have been carried out recently by Eiguren et 
al. [H,[35| for (i) the Shockley surface-state hole (n = 0) 
at the f point of Al(100) and the (111) surfaces of the no- 
ble metals Cu, Ag, and Au (34[, and (ii) the first (n = 1) 
image-state electron at the f point of the (100) sur- 
faces of Cu and Ag jUJ ; these calculations are based on 
the use of (i) Thomas-Fermi screened Ashcroft electron- 
ion pseudopotentials, (ii) single-particle states obtained 
by solving a single-particle model one-dimensional (ID) 
Schrodinger equation, and (iii) a simple force-constant 
phonon model calculation that yields a phonon spectrum 
in good agreement with experimental data. 

A summary of the results reported by Eiguren et 
al. [34|, [HI is presented in Table fl] Electron-phonon 
linewidths are particularly relevant in the case of surface- 
state holes with energies very near the Fermi level, in 
which case the contribution from e-e interactions is very 
small. In the case of image states, whose energies lie typ- 
ically a few electronvolts above the Fermi level, the e-ph 
linewidth is found to be T ep < 1 meV, thereby showing 
the negligibly small role of phonons in the electron dy- 
namics of image-potential states. 



III. ELECTRON-ELECTRON COUPLING 

Let us consider an arbitrary many-electron system of 
density «o(r). I n the framework of many-body theory, 
the e-e linewidth (or decay rate) T ee of a quasiparticle 
(electron or hole) that has been added in the single- 
particle state ipi(r) of energy is obtained as the projec- 
tion of the imaginary part of the self-energy S(r,r';e.i) 
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over the quasiparticle-state itself [H, H3] : 

T ee = ^2 j dr ydrV*(r)ImS(r,r';ei)^-(rO, (4) 

where the =F sign in front of the integral should be taken 
to be minus or plus depending on whether the quasipar- 
ticle is an electron (e^ > ep) or a hole (e.; < £f), respec- 
tively, Ef being the Fermi energy. Alternatively Eq. ((4]) 
can be written as follows 



2 r E *+° 

T ee =T-j q+ / dr I dr'. 9 ° (r, r'; e) ImS(r, r'; Si ), 



where 



(5) 



g°(r, v';e)= i - {g°(t>, r; s) - [G°(r, r'; e )] *} , (6) 

G°(r, r';e) being the one-particle Green function of a 
nonintcracting many-electron system: 



e — Sf + vr\ 



(7) 



Here, "0/( r ) an d £f represent the complete set of eigen- 
functions and eigenvalues of a one-particle hamiltonian 
describing the noninteracting many-electron system. 



A. Self-energy: G°W and G°W° approximations 

To lowest order in a series-expansion of the self- 
energy in terms of the frequency-dependent time- 
ordered screened interaction W (r , r' ; ui) , the self-energy 
E(r,r';£i) is obtained by integrating the product of 
the interacting Green function G(r, r';Ei — ui) and the 
screened interaction W(r, r'; ui), and is therefore called 
the GW self-energy. If one further replaces the inter- 
acting Green function by its noninteracting counterpart 
G°(r,r';£i — ui), one finds the G°W self-energy. For the 
imaginary part, one can write 



ImX G ° w (r, r'; e 4 ) = )_J ' ^(r')Im^(r, r'; | Ei - e/ |) ^(r), 
/ 

where the prime in the summation indicates that the 
sum is extended, as in Eq. ([7]), over a complete set of 
single-particle states V/( r ) 01 energy ej but now with 
the restriction sp < £f < £% or £i < £f < £f- In 
terms of the one-particle noninteracting Green function 
G°(r, r';e), one finds 

ImE GOw/ (r,r'; £l ) = - /'"* de Im 5 °(r, r'; \e t - e\) 



(9) 



n Jo 

x lmW(r,r';e), 



where g°(r,r';e) is given by Eq. ([5]). Introducing either 
Eq. © or Eq. © into Eq. Q or Eq. ©, one finds an 



expression for the e-e linewidth that exactly coincides 
with the result one would obtain from the lowest-order 
probability per unit time for an excited electron or hole 
in an initial state ipi (r) of energy to be scattered into 
the state V'/CO of energy £/ by exciting a Fermi system of 
interacting electrons from its many-particle ground state 
to some many-particle excited state (38) . 

The screened interaction W(r,r';ui) entering Eqs. ([5]) 
and ([9]) can be rigorously expressed as follows 



W(r,r';w) = v(r, r') + J dri J dr 2 v(r,ri) 



x x(ri,r 2 ;w)u(r 2 ,r'), 



(10) 



v(r,r') representing the bare Coulomb interaction and 
x(r, r';u>) being the time-ordered density-response func- 
tion of the many-electron system, which for the positive 
frequencies (ui > 0) entering Eqs. ([5]) and ^ coincides 
with the retarded density-response function of linear- 
response theory. In the framework of time-dependent 
density- functional theory (TDDFT) 39], the exact re- 
tarded density-response function is obtained by solving 
the following integral equation [40l | : 

X(r,r';w) = x°(r,r';w) + J dn / dr 2 x°(r,ri;w) 

x {«(n , r 2 ) + r c [n ] (n , r 2 ; w)} x(r 2 , r'; u). (11) 

Here, x°( r , r '; w ) denotes the density-response function 
of noninteracting Kohn-Sham electrons, i.e., independent 
electrons moving in the effective Kohn-Sham potential of 
density- functional theory (DFT) [Il[ : 

X°(r,r» = ^ £(/;-/;) 

^(r)y,*(r)^(r'M*(r') 
x : j (12) 

UJ — Ej + Si + If] 

where represents a normalization volume, fi are Fermi- 
Dirac occupation factors [which at zero temperature take 
the form /j = Q(sf — £i), Q(x) being the Heaviside 
step function], and V^( r ) an d £, represent the eigenfunc- 
tions and eigenvalues of the Kohn-Sham Hamiltonian of 
DFT. The other ingredient that is needed in order to 
solve Eq. (fTTj) is the xc kernel f xc [rio](r,r';u)), which 
is the functional derivative of the unknown frequency- 
dependent xc potential 14 c [n](r, uj) of TDDFT, to be 
evaluated at no(r). 

In the random-phase approximation (RPA), 
f xc [no](r, r'; uj) is set equal to zero: 

X flpA (r,r'; W ) = x°(r,r'; W ) + j dr, | dr 2 x °(r, r i; uj) 

x W (r 1 ,r 2 ) X ^(r 2 ,r'; W ), (13) 
and the screened interaction W(r, r'; ui) is replaced by 

W°{r,r';u) = v(r,r') + J dr r J dr 2 v(r,ri) 



x X rpa {v u v 2 ;lo)v{v 2 ,v% (14) 

or, equivalently, 

W°(r,r';uj) = u(r,r') + J dr x / dr 2 «(r,n) 

x x°(r 1 ,r 2 ;w)W°(r 2 ,r';w), (15) 

which yields the so-called G°W° (or G°W-RPA) self- 
energy: 

ImS G0 ^(r,r'; £i ) = £ ' ^(r')Im W°(r> r'; | e< - £/ |) 
/ 

x ^(r), (16) 

or, equivalently: 



ImE 



G°W° 



i r\ e i~ £ f\ 

(r,r';£i) = - / de Img°(r, r'; |e l - e|) 



(17) 




FIG. 2: Probability amplitudes <j>n{z) and energies e n of 
image-potential-induced bound states (n = 1, 2, 3) outside an 
infinitely repulsive solid surface occupying the z < space 
(shaded area). The thick solid line represents the classical 
image potential of Eq. (f2Tj) . 



x lmW°(r,r';e). 



C. Surface-state wave functions 



1. Simple models 



B. Self-energy: GWT approach 

The xc kernel f xc [no](r, r'; u>) entering Eq. (fTTj) . which 
is absent in the RPA, accounts for the presence of an xc 
hole associated to all screening electrons in the Fermi sea. 
Hence, one might be tempted to conclude that the full 
G°W approximation [with the formally exact screened 
interaction W of Eq. (flQ|) ] should be a better approxima- 
tion than its G°W° counterpart (with the screened inter- 
action W evaluated in the RPA). However, the xc hole as- 
sociated to the excited electron (or hole) is still absent in 
the G°W approximation. Therefore, if one goes beyond 
RPA in the description of W, one should also go beyond 
the G°W approximation in the expansion of the electron 
self-energy in powers of W. By including xc effects both 
beyond RPA in the description of W and beyond G°W 
in the description of the self-energy [52, EH, [3] , the so- 
called GWT approach yields a self-energy that is of the 
G°W form: 

Im£ GH/r (r, r'; Ei ) = ^(r')lm W(r, r'; | £i - e/ |) ^(r), 



or, equivalently: 



ImE 



GWT 



(r,r';ei) = - 



[si- 



de Im<7°(r, r'; 



ImW>(r,r';e), 



-e\) 
(19) 



but with the actual screened interaction W(r,r';u>) en- 
tering Eq. (jSJ) being replaced by a new effective screened 
interaction 

W{t,T f ;u)=v(T,T f )+ / dri / dr 2 {«(r,ri) 



+ r c [n ](r,r i;w )} X (n, r 2 ; u) v(r 2 , r'), (20) 
which includes all powers in W beyond the G°W approx- 



a. Outside the solid. Image states are quantum 
states trapped in the long-range image-potential well out- 
side a solid surface that presents a band gap near the 
vacuum level. In the case of a metal that occupies the 
half-space z < 0, the asymptotic form of the potential 
experienced by an electron in the half space z > is the 
classical image potential 



v '< 2) " "5 



(21) 



If one assumes (i) translational invariance in the plane 
of the surface and (ii) that due to the presence of a 
wide band gap at z < the solid surface is infinitely 
repulsive, i.e., V(z) — > oo at z < (45|, then one easily 
finds that the solutions of the corresponding one-particle 
Schrodinger equation represent a Rydberg-like series of 
image-potential induced bound states (see Fig. [2]) of the 
form: 



V^„(r) = 4=0nWe lk - r ", 
VA 



(18) with energies 



£n + k 2 /2, 



where 



and 



(j) n (z) ~ Z(f>] 



hydrogen 



(z/4), n = l,2,. 



e n = -- 



32n 2 



71 = 1,2,. 



(22) 



(23) 



(24) 



(25) 



imation. 



^hydrogen^ representing the well-known wave functions 
of all possible s-like (I = 0) bound states of the hydrogen 
atom. Here, r 1 1 and k represent the position and the 
wavevector in the xy surface plane. 
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FIG. 3: Probability amplitudes <j> n (z) of (a) Shockley (n = 0) 
and (b) image (n = 1) states on Cu(lll). The dotted lines 
represent the results one obtains by matching at z = the 
wave function of Eqs. (|26[) and (|27|l to a mere exponential (for 
the Shockley state) or to the hydrogenic form of Eq. (I24|l with 
n — 1 (for the n = 1 image state). The solid lines represent 
the results one obtains by solving a ID Schrodinger equation 
with the model potential of Chulkov et al. [49l |. 

b. Inside the solid. In the interior of the solid (z < 
0), both image and Shockley surface states can be de- 
scribed within a two-band approximation to the nearly- 
free-electron (NFE) band structure of the solid [46[ . As- 
suming translational invariance in the plane of the surface 
and for a gap that is opened by potential Fourier com- 
ponents corresponding to reciprocal lattice vectors that 
are normal to the surface, surface-state wave functions 
within the crystal band gap take the form 

4>(z < 0) re e A cos(Gz + <5). (26) 

Here, G represents the limit of the Brillouin zone in the 
direction normal to the surface, and 

where V g and e denote the energy gap and the surface- 
state energy with respect to the mid gap, respectively, 
and S represents a phase shift which in the presence of a 
Shockley- inverted band gap [47j varies from — tt/2 for a 
surface-state energy e at the bottom of the gap to for 
a surface-state energy at the top of the gap. Matching 
at z — to a wave function of the hydrogenic-like form 
of Eq. (in the case of image states) or to a mere 
exponential (in the case of Shockley states) [48| , one finds 
the wave functions 4> n (z) plotted by dotted lines in Fig. [3] 
for Cu(lll). 

2. One-dimensional model potentials 

Still assuming translational invariance in two direc- 
tions, i.e., assuming that the charge density and one- 
electron potential are constant in the plane of the surface, 



Chulkov et al. 49] devised a simplified model that allows 
for realistic calculations while retaining at the same time 
the essential physics of electron and hole dynamics at 
solid surfaces. In the bulk region, this one-dimensional 
(ID) model potential is described by a cosine function 
which opens the energy gap on the surface of interest, 
the position and amplitude of this function being cho- 
sen to reproduce the energy gap observed experimen- 
tally and/or obtained from first-principles calculations 
at the r point. At the solid-vacuum interface, it is repre- 
sented by a smooth cosine-like function that reproduces 
the experimental energy of the Shockley surface state. 
Finally, in the vacuum region this ID potential merges 
into the long-range classical image potential of the form 
of Eq. (|21|) in such a way that the experimental binding 
energy of the first image state is reproduced. This model 
potential has been constructuted for several metal sur- 
faces [5(| , and has been used widely for the investigation 
of electron and hole dynamics in a variety of situations. 

The n = and n = 1 eigenfunctions of a single- 
particle ID hamiltonian that includes the model poten- 
tial of Chulkov et al. for Cu(lll) are plotted in Fig. [3] 
by solid lines, together with the NFE Shockley (n = 0) 
and first image-state (n = 1) wave functions described 
in the preceding section. In the bulk region, these wave 
functions coincide with the approximate NFE wave func- 
tions (represented in Fig. [3]by dotted lines); however, in 
the vacuum region the n = 1 hydrogenic-like wave func- 
tion of Eq. (j2~4"|) appears to be too little localized near 
the surface. The n = 1 eigenfunction of Chulkov's ID 
hamiltonian for Cu(100) [51[ was found to reproduce ac- 
curately the average probability density derived for that 
image state by Hulbert et al. [52| from a first-principles 
calculation. 

The assumption that the charge density and one- 
electron potential are constant in the plane of the surface 
is valid for the description of image states, since their 
wave functions lie mainly at the vacuum side of the sur- 
face and the electrons move, therefore, in a region with 
little potential variation parallel to the surface. Shockley 
and bulk states, however, do suffer a significant poten- 
tial variation in the plane of the surface. In order to 
account approximately for this variation, the original ID 
model potential of Chulkov et al. [HI, which had been 
introduced to describe the projected band structure at 
the r point, was modified along with the introduction 
in Eq. (j23|) of a realistic effective mass for the disper- 
sion curve of both bulk and surface states [H[. Within 
this model, however, all Shockley states have the same 
effective mass, so the projected band structure is still in- 
accurate, especially at energies above the Fermi level, as 
shown in Fig. H for Cu(lll). 

As an alternative to the ID model potential of Chulkov 
et al. (4{| , Vergniory et al. [H[ introduced a k-dependent 
ID potential that is constructed to reproduce the actual 
bulk energy bands and surface-state energy dispersion 
obtained from 3D first-principles calculations, thereby al- 
lowing for a realistic description of the electronic orbitals 




0.2 0.3 

k (a.u.) 



FIG. 4: Dispersion of the Cu(lll) Shockley surface state 
(thick solid line), as obtained from 3D ab initio calcula- 
tions 53]. Shaded areas represent areas outside the band gap, 
where bulk states exist. The thin solid and dotted lines repre- 
sent approximate energy dispersions of the Shockley surface 
state and the bottom of the projected band gap, respectively, 
as obtained from an equation of the form of Eq. (I23p with the 
Cu(lll) f-point Shockley-state energy (with respect to the 
Fermi level) e = — 0.44eV and the effective mass m = 0.42 
(thin solid line), and with the Cu(lll) F-point bottom-of-the- 
gap energy (with respect to the Fermi level) e = —1.09 eV and 
the effective mass m — 0.22 (thin dotted line). 



beyond the T point: 

Uk + 2Vk cos(2w z/a s ), z < z k 



V k {z) 



(28) 



z > z k . 



Here, U k and V k are fitted to the bulk energy bands, a s 
represents the interlayer spacing, $ is the experimentally 
determined work function, and the matching plane z k is 
chosen to give the correct surface-state dispersion. 

The abrupt ID step model potential of Eq. (f2"8")) . which 
does not account for the image tail outside the surface, 
could not possibly be used to describe image states. How- 
ever, it has proved to be accurate for the description of 
Shockley surface states, which are known to be rather in- 
sensitive to the actual shape of the potential far outside 
the surface; indeed, the model potential of Eq. ([28| is 
found to yield a surface-state probability density |</>(z)| 2 
at the band edge (f point, i.e., k = 0) of the Shockley 
surface-state band of Cu(lll) that is in reasonably good 
agreement with the more realistic surface-state probabil- 
ity density obtained at f from the ID model potential of 
Chulkov et al. [49j , as shown in Fig. [5] Both probability 
densities coincide within the bulk, although the proba- 
bility density obtained from the step model potential of 
Eq. (|28[) appears to be slightly more localized near the 
surface, as expected. For the overlap integral one finds 
< <f>i\(/>2 >= 0.99, 0i and (j>2 being the Shockley proba- 
bility amplitudes leading to the probability densities rep- 
resented in Fig. [5] by solid and dashed lines, respectively. 




z (a.u.) 

FIG. 5: Probability-density |0(z)| 2 of the Shockley surface 
state at the center of the surface Brillouin zone (F point) of 
Cu(lll), as obtained with the use of two different ID model 
potentials: (i) the model potential of Chulkov et al. [4$J, which 
includes the image tail outside the surface but fails to repro- 
duce the actual band structure beyond the F point (solid line), 
and (ii) the model potential of Eq. (|28J1 . which does not in- 
clude the image tail outside the surface but is constructed to 
reproduce the actual bulk energy bands and surface-state en- 
ergy dispersion beyond the F point (dashed line). Full circles 
represent the atomic positions of Cu in the (111) direction. 
The geometrical (jellium) electronic edge [z = 0) has been 
chosen to be located half an interlayer spacing beyond the 
last atomic layer. 



D. Screened interaction 

The retarded counterpart of the density-response func- 
tion entering Eq. (fTT)|) , which in the framework of TDDFT 
can be obtained rigorously by solving the intergral 
Eq. (fTT|) . yields, within linear-response theory, the elec- 
tron density Sn(r, ui) induced in a many-electron system 
by a frequency-dependent external potential (jf xt {r,uj): 



5n{r,w) = J dr'x(r, r' ; tv) 4> ext (r' , uj) . 



(29) 



Hence, the retarded counterpart of the screened interac- 
tion W(t,t';oj) of Eq. (fTU)) yields, within linear-response 
theory, the total potential <j>(r,u>) of a unit test charge 
at point r in the presence of an external test charge of 
density n ext (r, u>); 

(f>(r, w) = J dr' W(r, r'; u) n ext (v', w), (30) 

which can also be expressed as follows 

0(r,w) = J dr' 'e-^r/jw) (f> ext (r, u), (31) 

with 

e-^r, r'; u) = 6(r -r') + J dr" «(r - r") X (r", r'; to). 

(32) 
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This is the so-called inverse dielectric function of the 
many-electron system, whose poles dictate the occur- 
rence of collective electronic excitations. 

1. Classical model 

In a classical model consisting of a scmiinfinite solid at 
z < characterized by a local (frequency-dependent) di- 
electric function e(u>) separated by a planar surface from 
a semiinfinite vacuum at z > 0, the total potential </>(r, u>) 
at each medium is a solution of Poisson's equation 

V 2 0(r,c)= n e * 4 (r, W ), (33) 

ti being e(u>) or 1 depending on whether the point r is 
located in the solid or in the vacuum, respectively. Hence, 



the screened interaction W(r,r';u>) entering Eq. (f30]) is 
a solution of the following equation: 



V 2 W(r,r» = — 5(r-r'). (34) 

Imposing boundary conditions of continuity of the po- 
tential and the normal component of the displacement 
vector at the interface, one finds 



W(r,r';uj)= dqe^^-^ W(z,z?;q,w), (35) 



where 



W{z,z';q,u) 



2tt 

q 



e -q\z^z'\ + g ^ e -q(\ Z \+z'\) / e (u) t Z<0,Z'<0, 

2 g(w) e - qz - z '\ I [e(w) - 1] , z< < 0, z> > 0, 

e — _ g (u) e -9(W+l 2 'D, z > 0, z' > 0, 



(36) 



z < (z > ) being the smallest (largest) of z and z 1 
being the classical surface-response function: 



e(w) + 1 ' 



and g 



(37) 



2. Specular-reflection model (SRM) 



An inspection of Eqs. (|3"B"|) and (|3"T|) shows that the 
screened interaction W(z, z';q,cu) has poles at the clas- 
sical bulk- and surface-plasmon conditions dictated by 
e(uj) = and by e(u}) + 1 = 0, respectively [5ij . Since e-e 
inelastic linewidths of Shockley and image states are typ- 
ically dominated by the excitation of electron-hole (e-h) 
pairs and not by the excitation of plasmons (whose en- 
ergies are typically too large) (5f| , the classical screened 
interaction of Eq. (f36|) (which obviously does not account 
for the excitation of e-h pairs) is of no use in this context. 



A simple scheme that gives account of the excitation of 
e-h pairs, and has the virtue of expressing the screened in- 
teraction W(z, z'\ q, lj) in terms of the dielectric function 
e(q,u>) of a homogeneous electron gas representing the 
bulk material, is the so-called specular-reflection model 
reported independently by Wagner and by Ritchie 
and Marusak (58|. In this model, the semi- infinite solid 
is described by an electron gas in which all electrons 
are considered to be specularly reflected at the surface, 
thereby the electron density vanishing outside. One finds: 



W(z,z';q,u) 



e s (z - z') + e s {z + z 1 ) -2g(q, u) e ^ € ^'\ z < , z> < 0, 



2tt 
— < 

q 



I \ e s(z < ) -„ z > 

m,u) 1 _ o e q , 

_ e -q\z-z'\ - 5 (g, W )e-«^+ z '), 
I 



z< <0,z> > 0, 
z > 0, z' > 0, 



(38) 



where the surface response function is now given by the 
following expression: 

l-e s (q,u) 



with 



e s (z;q,uj) = 



l + e2(«,o;)' 



(39) 



(40) 
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e°(<2S^) = e s (z = 0;q,u), 



(41) 



and Q = y/q 2 + g|. If the Q-dependence of the actual 
e(Q,w) dielectric function of a homogeneous electron gas 
is ignored, the SRM screened interaction of Eq. ([55)) re- 
duces to the classical screened interaction of Eq. (1361) . 

The inverse dielectric function e~ 1 (Q,uj) entering 
Eq. (|40|) represents the 3D Fourier transform of the in- 
verse dielectric function e _1 (r,r',cj) of a homogeneous 
electron gas. From Eq. (|3"2")l , one finds: 



(42) 



-\Q,u) = l + v QX (Q,u), 



where x(Q,w) represents the 3D Fourier transform of 
the density-response function x( r , r '; w ), and vq is the 
3D Fourier transform of the bare Coulomb interaction 
v(r,r'): vq = 4tt/Q 2 . 

In the framework of TDDFT, one uses Eq. (TTTj) to find 



x{Q,oj) = x\Q,u J ) + x {Q^) 

x{v Q + f xc (n ;Q,uj)} x(Q,w), 



(43) 



with x°(Qi UJ ) and f xc ( n o]Qj ljJ ) being the 3D Fourier 
transforms of the nonintcracting density-response func- 
tion x°( r - r '; w ) and the xc kernel f xc [no)(r, r'; uj) of a 
homogeneous electron gas of density uq. For a homo- 
geneous electron gas, the eigenfunctions ipi(r) entering 
Eq. (fT2|) are all plane waves; thus, the integrations can 
be carried out analytically to yield the well-known Lind- 
hard function x°(<3, w ) [321 . If one sets the xc kernel 
f xc (no;Q,uj) equal to zero, the introduction of Eq. 
into Eq. (|42[) yields the RPA dielectric function 



< rpa (Q,uj) = 1-vqx°(Q,uj), 



(44) 



which is easy to evaluate. 



3. ID self-consistent scheme 

For an accurate quantal description of the electronic 
excitations that occur in a semi-infinite solid, we need to 
consider the true self-consistent density-response func- 
tion x(r, r';ui) entering Eqs. (fTTTj) and (12"01) . 

Assuming translational invariance in the plane of 
the surface, one can still define the 2D Fourier trans- 
forms W(z, z'\q,uj) and W(z, z';q,uj), which according 
to Eqs. (jlOp and (|20|) can be obtained as follows 



W(z, z'; q,uj) = v(z, z'; q) + J dz\ J dz 2 v(z, Z\;q) 
x x( z i, z 2;q,u)v(z 2 ,z';q), (45) 



and 



W(z, z'\ q, w) = v(z, z'\ q) + J dzi J dz 2 {v(z,z%;q) 
+ f xc [no}(z,z 1 ;q,u)} x{zi,z 2 ;q,uj)v{z 2 ,z',q), (46) 



where v(z, z'; q) is the 2D Fourier transform of the bare 
Coulomb interaction u(r,r'): 



v(z,z';q) = — e -i\z-z'\, 

q 



(47) 



f xc [n](z, z'\ q,u>) is the 2D Fourier transform of the xc 
kernel f xc [n](r, r'; uj), and x( z > z '\ 1i w ) denotes the 2D 
Fourier transform of the interacting density-response 
function x(r,r';o;). In the framework of TDDFT, one 
uses Eq. (fTTj) to find: 



X(z,z';q,uj) = x°(z,z';q,uj) + J dz x J dz 2 x°(z,zx;q,uj) 

x {v(z 1 ,z 2 ;q) + f xc [n ]( Zl , z 2 ;q,uj)} x( z 2, z ';q,u), (48) 

where x°( z ? z '\ 1i w ) denotes the 2D Fourier transform of 
the noninteracting density-response function x (r, r'jw). 
Using Eq. (fl"2"|) , and noting that the single-particle or- 
bitals ipi(r) now take the form 



^k,i(r) = < (z)e ut -'n 



(49) 



one finds: 



X°{z,z';q,uj) = — ^2 M z )(t>*3( z )<t>j( z ')(t>U z ') 



E 



/k,i — /k+q,j 



Ek,i - £k+q,j + w + iv ' 



(50) 



where /k,i are Fermi-Dirac occupation factors [which at 
zero temperature take the form /k,i — 0(sf — Ek,i)], and 



k 2 

E\t t i — Si + — , 



(51) 



the single-particle orbitals 4>i(z) and energies £j now be- 
ing the eigenfunctions and eigenenergies of a ID Kohn- 
Sham hamiltonian. In order to account for the actual 
band structure of sp electrons near the surface of sim- 
ple and noble metals, <pi{z) and have been succes- 
fully taken to be the solutions of the ID single-particle 
Schrodinger equation of Chulkov et al. [49[ described in 
the previous section. 



4- Asymptotics 

For z and z' coordinates that are far from the surface 
into the vacuum, where the electron density vanishes, 
Eq. J45]) takes the form [6(J 

W(z, z'; q, u) = v(z, z'; q) - — g(q, uj), (52) 

with the surface-response function g(q,u>) now being 
given by the general expression [6lT | : 

g(q,uj) = -— J dzx J dz 2 e^ +z ^ X (zi,z 2 ; q,u). 

(53) 
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Persson and Anderson [62j and Persson and 
Zaremba |63| investigated the structure of the so- 
called surface loss function Img(q,uj) for small q and uj. 
Persson and Zaremba found the following approximate 
result (6|: 



and m 



surf 



(Imp) 



bulk 



(Img)i 



(54) 



where (Iuig) sur f and (lmg)buik represent contributions 
from surface and bulk excitation of c-h pairs, and 
(Im g)int represents the contribution to the surface loss 
function coming from the interference between the bulk 
and surface excitations: 



and 



q uj 

(lmg) surf = 2£- , 

Kf UJ p 



{Img)buik = ^mlpt (j^- ) ifG[q). 



8 m opt ( u\ h 
{img) int = j T — — % , 2 - 



(55) 



(56) 



(57) 



with rj = Lu/(qkp) and 



for rj < 1, 



G(rj) 



l-(l + ir7- 2 )(l-7 ? - 2 ) 1 / 2 , for?y>l. 



(58) 

Here, k F and uj p represent the Fermi momentum and 
the plasmon frequency, respectively: kp — 37r 2 no and 
ujp = 47rno, uq being the mean electron density. The 
values of £, m opt , and h are given in Ref. (63j . The surface 
contribution of Eq. (|55p had already been reported in 
Ref. [62j, the bulk contribution of Eq. (p)6"|) differs from 
that used in Ref. [62j by the factor of the optical mass 
and a factor of 5 which had been missed previously, and 
the contribution (lmg) int had been neglected in Ref. [6^ ]. 



5. The xc kernel f xc [n](z, z';q,u>) 

Several approximations can be used to evaluate the 
unknown xc kernel f xc [no](z, z'\ q, u>) entering Eqs. J46|) 



a. Random-phase approximation (RPA). Nowadays 
one usually refers to the RPA as the result of sim- 
ply setting the xc kernel f xc [n](z, z'\ q, uf) equal to zero: 
f xc [n] (z, z'; q, w) = 0, but still using in Eqs. |@5J and f5T]l 
single-particle Kohn-Sham states and energies 4>i(z) and 
Si that go beyond the Hartree approximation. 

b. Adiabatic local-density approximation (ALDA). 
If one assumes that dynamic electron-density fluctua- 
tions are slowly varying in all directions, the xc kernel 
f xc [n] (z, z'; q, uj) is easily found to be given by the fol- 
lowing expression [6l| : 



f xc [n}(z,z';q,u;) = f xc (n = n(z);Q = 0,lu = O)S(z-z'). 

(59) 

Here, f xc (n = n{z); Q,uj) is the 3D Fourier transform of 
the xc kernel of a homogeneous electron gas of density 
n equal to the local density n(z), which in the limit as 
Q — > and u> — > is known to be the second deriva- 
tive of the xc energy per particle e xc (n). One typi- 
cally uses parametrizations [64] of the diffusion Monte 
Carlo (DMC) xc energy e xc reported by Ceperley and 
Alder [6|. 

c. Adiabatic nonlocal approximation (ANLDA). 
The investigation of short-range xc effects in solids has 
focused to a great extent onto the homogeneous electron 
gas [gfjj]. Hence, assuming that the unperturbed density 
variation [n(z) — n(z')] is small within the short range 
of f xc [n](z, z'; q,u>), one can adopt the following average 
adiabatic nonlocal approximation [67t l68l . |69] ] : 



f xc [n](z, z'; q, lo) = f xc ([n(z) + n(z')\ /2; z, z'; q,LJ = Q), 

(60) 

where f xc (n; z, z'\ q,w) represents the ID Fourier trans- 
form of the xc kernel f xc (n; Q, u>) of a homogeneous elec- 
tron gas of density n. A parametrization of the accu- 
rate DMC calculations reported by Moroni et al. [7(| 
for the static (u = 0) Q-dependent nonlocal xc kernel 
f xc {n; Q, uj — 0) that satisfies the well-known small and 
large-wavelength asymptotic behaviour was carried out 
by Corradini et al. (CDOP) 71]. An explicit expression 
for the 2D Fourier transform of the CDOP parametriza- 
tion of f xc (n; Q, uj = 0) was reported in Ref. [69l |: 



f xc (n;z,z';k) = 



4ire 2 C 
hp 



S(z) 



2ne 2 B 
V9 k F + fc2 



Vfffe|+fe 2 |5| 



7,3 



20 -kj 
4/3 2 



-hp 



-P[klz 2 /4/3 2 + fe 2 jk\ 



where C, B, g, a, and (3 are dimensionless functions of 
the electron density (see Ref. [Tljj ) , k F = (S^n) 1 ^ , and 
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Calculations of the frequency dependence of the xc ker- 
nel of a homogeneous electron gas have been carried out 
mainly in the limit of long wavelengths [ll, 0, [ll, 
17(1 l77f . but work has also been done for finite wave vec- 
tors |78l. [79l . [80L [8l| . Approximate expressions for the 
frequency-dependent xc kernel of inhomogeneous systems 
have been reported in Refs. [U 



6. d— electron screening 

The ID self-consistent scheme described above has 
proved to be appropriate for the description of the 
screened interaction of sp electrons in simple and noble 
metals. It has been argued, however, in the past that a 
realistic first-principles description of the electronic band 
structure is of key importance in the determination of 
the inelastic lifetime of bulk electronic states in the no- 
ble metals, due to the participation of d electrons in the 
screening of e-e interactions |8a |. 

Following the scheme originally developed by Liebsch 
to describe the anomalous dispersion of surface plasmons 
in Ag [87} , Garcfa-Lekue et al. [56| accounted for the pres- 
ence of occupied d-bands in the noble metals by assuming 
that sp valence electrons are embedded in a polarizable 
background at z < zq characterized by a local dielectric 
function ed(u>)- Within this model, the bare Coulomb 
interaction v(z,z';q) entering Eq. l[38"|) is replaced by a 
modified (<i-screened) Coulomb interaction v'(z,z';q,u)) 
whose 2D Fourier transform yields [88| 



v'(z,z';q,u) 



2tt 



qe d (z',uj) 
x aa(u})e 



«l^ z l+sgn(z d -z') 

(62) 



q\ z - z d\ e -q\zd-z'\ 



where 



and 



ed(z,uj) 



€d(w), Z < Z d 



Z > Z d 



Od = 



edM + i" 



(63) 



(64) 



The first term in Eq. (f6"2"|) is simply the 2D Fourier trans- 
form of the bare Coulomb interaction [see Eq. (|47pj. but 
now screened by the polarization charges induced within 
the polarizable background. The second term stems from 
polarization charges at the boundary of the medium. 



x W StS ,(z,z';q,u>), 



(65) 



where q is a 2D wave vector in the surface Brillouin zone 
(SBZ), and g and g' denote 2D reciprocal-lattice vec- 
tors. According to Eq. (flQ|) . the 2D Fourier coefficients 
Wg. g '(z, z'; q, uj) are given by the following expression: 



W Stg '(z, z'\q,u>) =v g (z,z';q)6 gt g/ + J dz x J dz 2 

x v s [z, zi ; q) Xg, g ' (zi , z 2 ; q, w) V (z 2 , z'\ q), (66) 

where v s (z,z';q) denote the 2D Fourier coefficients of 
the bare Coulomb interaction v(r, r'): 



u g (z,z';q) = 



2tt 



q+g 



-|q+g||*-*'| 



(67) 



and Xg,g'( z i z'; q, ui) are the Fourier coefficients of the 
interacting density-response function x(r, r';cj). In the 
framework of TDDFT, one uses Eq. (fTT|) to find: 

Xg, g '(z,z';q,u) =x s , S '(z,z';q,Lj) + J dz x J dz 2 
xx gtf ,>{z,zi;q,w) x [v Sl (zi,z 2 ;q)S Sug2 
+/g J i C ,g 2 [ n o]( z i: z 2;q,w)] x g2 ,g'( z 2,z';q,w), (68) 

where x° s , s > (z, z'\ q, uj) and f™ g , [n ] (z, z'\ q, w) denote 
the Fourier coefficients of the noninteracting density- 
response function x°( r ! r '; ti; ) an d the xc kernel 
/ xc [no](r, r'; uj), respectively. Using Eq. (fT2"]) , one finds: 



SBZ 



/k,n /k+q,n' 



x(0 k ,„W|e-^ + s> r ll|0 k+q ^(^)) 
x(^ k+q ^(z')|e t(q+S ' ) ' r|l l^n(2) 



(69) 



the single-particle orbitals $>k,n( r ) and energies £k, n being 
the eigenfunctions and eigenvalues of a 3D Kohn-Sham 
Hamiltonian with an effective potential that is periodic 
in the plane of the surface. 



IV. RESULTS AND DISCUSSION 



A. Image states 



7. Periodic surface 

For a real periodic surface, one may introduce the fol- 
lowing Fourier expansion of the screened interaction: 

SBZ 

W{v,t'-uj) = i ^ ^ c >( q+g )- r iie- j(q+s 'Mi 
i g,g' 



The first quantitative evaluation of image-state life- 
times was reported in Ref. [891 ] . This calculation was 
carried out from Eqs. ((4]) and ([8]), with (i) the hydrogenic- 
like image-state wave function <f>i{z) of Eq. with 
n = 1 and no penetration into the solid, (ii) the bulk 
final state wave functions 4>f{ z ) obtained with the use of 
a step model potential, and (iii) two simplified models for 
the screened interaction: the SRM of Eq. (f3"8"|) with the 
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RPA for the bulk dielectric function, and the surface re- 
sponse function reported by Persson and Anderson [62| . 
In subsequent calculations the penetration of the image- 
state wave function into the crystal was allowed [90| , and 
the role that the unoccupied part of the narrow Shockley 
surface state on the (111) surfaces of Cu and Ni plays 
in the decay of the n = 1 image state on these sur- 
faces was investigated by Gao and Lundqvist [oi| . In 
this work, the image-state wave functions were also ap- 
proximated by hydrogenic-like wave functions of the form 
of Eq. (|24[) with no penetration into the solid, a simpli- 
fied parametrized form was used for the Shockley surface- 
state wave function, and screening effects were neglected 
altogether. A G°W° calculation of the imaginary part of 
the electron self-energy near a jellium surface was also re- 
ported [92[ , showing the key role that a full evaluation of 
this quantity may play in the description of surface-state 
lifetimes. 

The first self-consistent many-body calculations of 
image-state lifetimes on noble and simple metals were 
reported by Chulkov et al. [HI, [93| , and good agreement 
with the experimentally determined decay times [9J, l95l . 
l96j was found. In these calculations, all wave functions 
and energies were obtained by solving a single-particle 
Schrodinger equation with the physically motivated ID 
model potential of Ref. [49[ , and the electron self-energy 
was evaluated in the G°W° approximation. The poten- 
tial variation in the plane of the surface was considered 
later through the introduction of an effective mass [48j . 
and self-consistent calculations of the key role that the 
partially occupied Shockley surface state plays in the de- 
cay of image states on Cu(lll) were carried out [48| . 
The inclusion of xc effects was investigated in the frame- 
work of the GWT approximation, first with an adiabatic 
local-density description [97j and more recently with an 
adiabatic non-local description of the xc kernel [98| . 

The impact of xc effects on the imaginary part of the 
effective screened interaction of Eq. (|46|) in the vicinity of 
the (100) and (111) surfaces of Cu and Au is illustrated 
in Fig. [HI where ANLDA calculations of Im[W(^, z'\ q, ui)] 
(with full inclusion of xc effects) are compared to calcu- 
lations oflm[W(z,z';k,uj)] with (ANLDA) and without 
(RPA) xc effects. Exchange-correlation effects included 
in the effective screened interaction have two sources. 
First, there is the reduction of the screening due to the 
presence of an xc hole associated with all electrons in the 
Fermi sea [see Eq. (|48p], which is included in the calcu- 
lations represented in Fig. [5] by thick solid lines and also 
in the calculations represented by dotted lines. Secondly, 
there is the reduction of the effective screened interac- 
tion itself due to the xc hole associated with the excited 
electron or hole [see Eq. (|46p] . which is only included in 
the calculations represented in Fig.[S]by thick solid lines. 
These contributions have opposite signs, thereby bring- 
ing the ANLDA Im[W(z, z'\ k, lj)] (thick solid lines) back 
to the RPA Im[W(z, z'; k, w)] (thin solid lines). 

Figure exhibits G°W° : G°W, and GWT calcula- 
tions of the imaginary part of the k-resolved n = 1 



Cu(lOO) 0>=leV aj=2eV o>=3eV to=4eV 




-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 -10 -5 5 10 

z=z' (a.u.) 



FIG. 6: Imaginary part of the screened interaction 
W(z, z' ; q,ui) and the effective screened interaction 
W(z,z ;q,u), as a function of z = z for a fixed value 
of k (q = 0.5A ) and various values of u (w = l,2,3,4eV), 
in the vicinity of the (100) and (111) surfaces of Cu and Au. 
ANLDA calculations of Im [W"(^, z'; q, a>)l are represented 
by thick solid lines. RPA and ANLDA calculations of 
lm[W(z,z';q,Lj)] are represented by thin solid and dotted 
lines, respectively. The RPA lm[W(z, z' ; q,ui)] and the 
ANLDA Tai\W(z,z';q,u))] are nearly indistinguishable, as 
a result of a cancellation between the xc effects due to the 
presence of (i) an xc hole associated with all electrons in the 
Fermi sea and (ii) an xc hole associated with the excited 



Cu(lll)n=0 Cu(100)n=l 




z(a.u.) z(a.u.) 

FIG. 7: G°W° (G°VF-RPA), G°W, and GWT calculations 
of the imaginary part of the n = 1 image-state self-energy 
£(z,z';k = 0,£k), versus z, in the vicinity of the (100) sur- 
faces of Cu and Au. The solid circle represents the value 
of z in each case. GWT calculations (as obtained with the 
use of our ANLDA xc kernel) are represented by thick solid 
lines. G°W (also using our ANLDA xc kernel) and G°W° 
calculations are represented by dotted and thin solid lines, 
respectively. ALDA calculations, not plotted in this figure, 
are found to nearly concide with ANLDA calculations. 
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TABLE II: Reciprocal lifetimes, in linewidth units (meV), of 
the n = 1 image state (at T, i.e., k = 0) on Cu(100), tog ethe r 
with the most recent measurement reported in Ref. |l00l ]. 
All the single-particle wave functions and energies entering 
Eqs. @, ©, (|16|1 . and (I18|l have been chosen to be the eigen- 
functions and eigenvalues of the ID Hamiltonian of Chulkov et 
al. [49J] . Effective masses for all the single-particle energies en- 
tering Eqs. |[HJ|, (|16[) , and (JTSj have been set equal to either the 
free-electron mass (my = 1) or to realistic values (rrif 7^ 1). 
Various models have been considered for the description of 
the electron self-energy, as obtained from (i) Eq. (J8| with 
the screened interaction of Eq. ((38j (SRM), of Eqs. (l52"]) - ([58j) 
(PZ), and of Eqs. J45J and gSJ with f xc [n ]{z, z'\ q, w) set 
equal to zero (G°W°), (ii) Eq. (|16p with the screened in- 
teraction of Eqs. J45j and (l48l) and the ANLDA xc kernel 
f xc [n ](z, z'\ q, cu) (G°W), and (iii) Eq. Jl8j with the effective 
screened interaction of Eqs. (|46p and (|48[l and the ANLDA xc 
kernel f xc [n ]{z, z'; q, uj) (GWT). 

mf Self-energy Bulk Vacuum Mix Total Exp. 



=1 


SRM 


18 


3 


-4 


17 


=1 


PZ 




55 




55 


=1 


G°W° 


24 


14 


-16 


22 


/I 


G°W° 


7 


11.5 


-1 


17.5 


/I 


G°W 








24.5 


/I 


GWT 


6.5 


11.5 


-1 


17 



16 



TABLE III: As in Table [TTJ but for Cu(lll) and together 
with the reciprocal lifetime experimentally determined for this 
surface at low temperature, T = 25 K (951 ]. 

nif Self-energy Bulk Vacuum Mix Total Exp. 



=1 


SRM 


46 


12 


-22 


36 


=1 


PZ 




57 




57 


=1 


G°W° 


44 


47 


-54 


37 


/I 


G°W° 


32 


34 


-37 


29 


/I 


G°W 








43 


/I 


GWT 


30 


35 


-38 


28 



30 



image-state self-energy E(z, z';k = 0,£k), versus z, in 
the vicinity of the (100) surfaces of Cu and Au, with use 
(in the case of the G°W and GWT approximations) of 
the adiabatic nonlocal xc kernel (ANLDA) described in 
Sec. IIII D 51 This figure shows that as occurs in the case 
of the screened interaction xc effects partially compensate 
each other, leading to an overall effect of no more than 
5% percent. We note that, as anticipated in Ref. [98( for 
the case of Cu(lll), although the ALDA leads to spuri- 
ous results for the screened interaction our more realistic 
ANLDA kernel yields self-energies that esentially coin- 
cide with those obtained in the ALDA. 

Tables |TT] and IIIII exhibit the linewidth of the n — 1 
image state (at f , i.e., k = 0) on the (100) and (111) sur- 
faces of Cu, as obtained from Eq. ([4]) with (i) the image- 
state wave function of Chulkov et al. [49| described in 
Sec. IIII C 21 and (ii) various approximations for the self- 



energy: SRM, PZ, G°W°, G°W, and GWT. Contribu- 
tions to the linewidth are separated as follows 

Tee — ^bulk "T" T vac -{- T m i x , (^0) 

where T bu i k , T 

vaci and T m i x represent bulk, vacuum, and 
mixed contributions, respectively, as obtained by confin- 
ing the integrals in Eqs. (j4j to either bulk (z < 0, z' < 0), 
vacuum (z > 0,z' > 0), or mixed (z > 0, z' < or 
z < 0, z' > 0) coordinates. 

First of all, we set all effective masses equal to the 
free-electron mass, and focus on the role that an accurate 
description of the screened interaction plays in the cou- 
pling of image states with the solid, by comparing to the 
results obtained (within the G°W° approximation) with 
the use of the SRM screened interaction and (for the vac- 
uum linewidth) the screened interaction of Persson and 
Zaremba (PZ). We note that simplified jellium models for 
the evaluation of the screened interaction yield unrealistic 
results for the image-state lifetime. Bulk contributions to 
the linewidth are approximately well described within the 
SRM, small differences resulting from an approximate de- 
scription, within this model, of the so-called beqrenzung 
or boundary- effect first described by Ritchie [991 ]. How- 
ever, as quantum-mechanical details of the surface are 
ignored within this model, it fails to describe both vac- 
uum and mixed contributions to the decay rate. These 
quantum-mechanical details of the surface are approxi- 
mately taken into account within the PZ approach, but 
the PZ model cannot account for the coupling of the im- 
age state with the crystal that occurs through the pen- 
etration of the image-state wave function into the solid. 
Discrepancies between vacuum contributions obtained in 
this model and in the more realistic full G°W° approach 
appear as a result of the PZ model being accurate only 
for small q wave vectors and to frequencies. 

Now we account for the variation of the potential in 
the plane of the surface through the introduction of a re- 
alistic effective mass for all surface and bulk states. The 
effective masses of the n = 1 image state on Cu(100) and 
Cu(lll) are close to the free-electron mass (rrii = 1). 
Nevertheless, the effective mass of the n = Shockley 
surface state of Cu(lll) and the unoccupied bulk states 
in Cu(lll) and Cu(100), which all contribute to the de- 
cay of the n = 1 image state, considerably deviate from 
the free-electron mass; Tables |TT] and IIIII show that the 
impact of this deviation on the n — 1 image-state lifetime 
is not neglegible. 

As for the impact of short-range xc effects, which are 
fully incorporated in the framework of the GWT approx- 
imation, Tables [Til an d IIIII show that the overall impact 
of these effects is small and GWT reciprocal lifetimes are 
close to their G°W° counterparts. 

ALDA and ANLDA GWT calculations of the recipro- 
cal lifetimes of the n = 1 image state on Au(100), never 
reported before, are exhibited in Table [IV] For compari- 
son, this Table also shows G°W° and G°W calculations, 
which in the case of the G°W have been obtained by con- 
sidering (as within the GWT approximation) both local 
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TABLE IV: G°W°, G°W, and GWT reciprocal lifetimes, in 
linewidth units (meV), of the n — 1 image state (at T, i.e., 
k = 0) on Au(100). In the case of the G°W and GWT re- 
ciprocal lifetimes, both ALDA and ANLDA xc kernels have 
been considered. 



xc kernel 


G°W° 


G°W 


GWT 




30 






ALDA 




42.5 


31 


ANLDA 




43 


31 



TABLE V: Linewidth (inverse lifetime) of image states, in 
meV. The lifetime in fs (lfs = 10 _lj s) is obtained by noting 
that h — 658meVfs. The numbers in brackets represent the 
corresponding references. Electron-phonon linewidths are not 
included, since in the case of image states they are all expected 
to be below 1 meV [35|. 





G°W° 


GWT 


Exp. 


Li(110) 


37 [97] 






Cu(100) 


17.5 [97] 


17 [98] 


16.5 ±2.5 [94, 96] 








16 [100] 


Cu(lll) 


29 [97] 


28 [98] 


30 [95j a 








29 ±6 [1011'' 


Ag(100) 


12 [56] 




12 ± 1 [96] 


Ag(lll) 


36 [56] 




21 ±9/5 [102] 


Au(100) 


30 [103] 


31 c 




Pd(lll) 


30 [104] 




26 ±5/3 [1041 


Pt(lll) 


23 [105] 




25 ± 10/5 [105J 


Ni(100) 


33 [103] 




41 ± 19/10 [106] 


Ni(lll) 


44 [103] 




94 ±71/28 [107] 


Ru(1000) 


47 [1081 




59 [1091 


a At T = 25 K. 
6 At T = 100 K. 
c This work. 



(ALDA) and nonlocal (ANLDA) xc kernels. As occurs 
in the case of Cu [§8|, the results shown in Table [IV] indi- 
cate that (i) a realistic adiabatic nonlocal description of 
xc effects yields reciprocal lifetimes of image states that 
esentially coincide with those obtained in the ALDA, and 
(ii) the overall effect of short-range exchange and corre- 
lation is small, thereby GWT reciprocal lifetimes being 
close to their G°W° counterparts. 

In Table [Vj we compare self-consistent G°W° and 
GWT (as obtained with the ANLDA xc kernel) calcu- 
lations (with full inclusion of realistic values of the effec- 
tive mass of all bulk and surface states) with the exist- 
ing TR-2PPE data for the n = 1 image state at the f 
point (k = 0) on various simple, noble, and transition 
single-crystal surfaces. This table shows that G°W° and 
GWT calculations are both in good agreement with TR- 
2PPE measurements except in the case of the (111) sur- 
faces of the noble metal Ag and the transition metal Ni. 
The largest disagreement occurs in the case of Ni(lll), 
where the G°W° linewidth is smaller than the measured 
linewidth by approximately a factor of 2. This can be at- 
tributed to the necessity of a full ab initio description of 



the dynamical response of both sp and d electrons along 
the lines of Sec. HIT D 71 

The role that occupied d bands play in the dynam- 
ics of i mag e states on silver surfaces was investigated 
in Ref. [56J along the lines described in Sec. HHP 61 It 
was concluded that d electrons do play an important role 
as a consequence of the reduction (in the presence of d 
electrons) of the surface-plasmon energy that allows the 
opening of a new decay channel. No surface-plasmon de- 
cay channel is opened, however, in the case of the other 
noble-metal surafces (Cu and Au), since in the presence 
of d electrons the Cu and Au surface-plasmon energy is 
still too large. 



B. Shockley states 

G°W° calculations of the e-e inelastic lifetimes of ex- 
cited holes at the surface-state band edge of the (111) 
surfaces of the noble met als C u, Ag, and Au were first 
reported in Ref s. [29T| and [l08j ] within the ID scheme of 
Chulkov et al. [49| (see Sec. IIIIC 2|) . accounting for the 
potential variation in the plane of the surface through 
the introduction of a realistic effective mass for the dis- 
persion curve of both bulk and surface states. Within 
this model, however, all Shockley states have the same 
effective mass and the projected band structure is still 
inaccurate, especially at energies above the Fermi level, 
as shown in Fig. 2] for Cu(lll). As an alternative to 
the ID model potential of Chulkov et al. ;49|, Vergniory 
et al. |53J introduced the k-dependent ID potential of 
Eq. (|28[) that is constructed to reproduce the bulk en- 
ergy bands and surface-state energy dispersion obtained 
from 3D first-principles calculations. 

Table I VII shows a comparison between the G°W° cal- 
culations reported in Refs. [29|, Il08j ] and [H[ for the in- 
elastic lifetime of an excited hole at the band edge of 
the Shockley surface-state band of Cu(lll), as obtained 
from Eqs. (|4]) and ([5]) with the use of the ID scheme of 
Chulkov et al. [49[ and with the k-dependent ID model 
potential of Eq. (|2"8|) , respectively. The difference be- 
tween the surface-state lifeti me broadening of 25 eV re- 
ported in Refs. [2!| and jl08| ] and the more accurate life- 
time broadening of 19 eV reported in Ref. 15311 is entirely 
due to a more accurate description in Ref. [53j of (i) the 
projected band structure and (ii) the wave- vector depen- 
dence of the surface-state wave functions entering the 
evaluation of the self-energy. G°W and GWT calcula- 
tions were reported in Ref. [§8j], showing that as in the 
case of image states GWT linewidths are only slightly 
lower than their G°W° counterparts. 

At this point, we note that the linewidths of the 
Cu(lll) Shockley state at T based on the use of the two 
ID models of Sec. IIII C 21 to describe the initial surface- 
state wave function (at f) agree within less than 1 meV. 
The linewidths also agree within less than 1 meV when 
the actual surface-state dispersion (thick solid line of 
Fig. 2]) is replaced by the parabolic surface-state disper- 
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TABLE VI: G°W°, G°W, and GWF decay rates, in linewidth 
units (meV) , of an excited hole at the band edge of the Shock- 
ley surface-state band of Cu(lll) (Ei = —0.44 eV and ki = 0). 
The G°W° calculations have been performed either with the 
use of the k-dependent ID model potential of Eq. (|28[) . as 
reported in Re f . [ 5311 . or with the use of the ID sche me of 
Chulkov et al. [ifTas reported in Refs. [U and [i"0g| . The 
G°W and GWF calculations have been per formed with the 
use of the ID scheme of Chulkov et al. [40], as reported in 
Ref. [§||. The experimental linewidth has been taken from 
the STM measurements reported in Ref. [29l ] . The total decay 
rate r to tai includes the e-ph decay rate of 7 meV reported in 
Ref. [34|. Tinter and Lmtra represent interband and intraband 
contributions to the e-e decay rate r e - e ; these contributions 
come from the decay of the excited hole trhough the coupling 
with bulk states (interband contribution) or thorugh the cou- 
pling, within the surface-state band itself, with surface states 
of different wave vector k parallel to the surface (intraband 
contribution) . 





Tintcr 


Tintra 


Tc — c 


Ttotal 


G l 'W u [53] 


10 


9 


19 


26 


G°W° [29. 1081 


6 


19 


25 


32 


G°W [98] 






30.5 


37.5 


GWF [98] 






24.5 


31.5 


Experiment 








24 



TABLE VIII: Linewidth (inverse lifetime) of Shockley states, 
in meV. The lifetime in fs (lfs = 10 _15 s) is obtained by not- 
ing that h = 658meVfs. The numbers in brackets represent 
the corresponding references. The calculated values (r ca ; c ) 
are decomposed into e-e (r e _ e ) and e-ph (r e _ p h) contribu- 
tions. Since GWF e-e linewidths are found to be very close 
to their G°W° counterparts, only G°W° calculations of the 
e-e linewidth are included here. In the case of Be(0001), cal- 
culations have been performed either with the u se of the ID 
scheme of Chulkov et al. [4j|, as reported in Ref. [l 10( | . or via 
a fully ab initio scheme along the lines of Sec. IHI D 71 as re- 
ported in Ref. [l 13f | . In the case of Cu(lll), calculations have 
been performed either with the use of the k-dependent ID 
model potential of Eq. (|28|l . as reported in Ref. [53I ]. or with 
the use of the I D sch eme of Chulkov et al. [4{J, as reported 
in Refs. [U and [TosT ]. 





r 

A e — e 






Exp. 


Al(lll) 


336 [110J 


36 [110J 


372 


~ 1500 [111] a 


Mg(0001) 


83 [110J 


25 [110] 


108 


~ 500 [nij a 


Be(0001) 


280 [110] 


80 [113] 


360 






265 [113] 


80 [113] 


345 


350 [113] 


Cu(lll) 


25 f29. 1081 


7 [34] 


32 






19 [53] ~ 


7 [34] 


26 


24 [29] 


Ag(lll) 


3 [29. 1081 


4 [34] 


7 


6.5 [29] 


Au(lll) 


29 [29. 1081 


4 [34] 


33 


18 [29] 


"At room temperature. 



TABLE VII: G°W°, G°W, and GWF reciprocal lifetimes, in 
linewidth units (meV), of the n — Shockley state (at F, 
i.e., k = 0) on Au(lll). In the case of the G°W and GWF 
reciprocal lifetimes, both ALDA and ANLDA xc kernels have 
been considered. 



xc kernel 


G°W° 


G°W 


GWF 




29 






ALDA 




39.5 


30 


ANLDA 




40 


30 



sion of the form dictated by the thin solid line of Fig. [4] 
However, if one replaces in the calculation of Ref. [53| the 
wave- vector dependent surface-state orbitals obtained by 
solving a ID Schrodinger equation with the potential 
of Eq. (I28|) by thei r less accurate counterparts used in 
Refs. |29l] and |l08j . the lifetime broadening is increased 
considerably (from 19 to 25 meV), showing the impor- 
tant role that the actual coupling between initial and 
final states plays in the surface-state decay mechanism. 

ALDA and ANLDA GWT calculations of the recipro- 
cal lifetimes of the Shockley surface state on Au(lll), 
never reported before, are exhibited in Table IVIII For 
comparison, this Table also shows G°W° and G°W calcu- 
lations, which in the case of the G°W have been obtained 
by considering (as within the GWT approximation) both 
local (ALDA) and nonlocal (ANLDA) xc kernels. As oc- 
curs in the case of Cu [§8|, the results shown in Table fVTTl 
indicate that (i) a realistic adiabatic nonlocal description 
of xc effects yields reciprocal lifetimes of Shockley states 
that esentially coincide with those obtained in the ALDA, 



and (ii) the overall effect of short-range exchange and cor- 
relation is small, thereby GWT reciprocal lifetimes being 
close to their G°W° counterparts. 

The calculated and experimental linewidths of Shock- 
ley states at the T point of a variety of simple and noble 
metal surfaces are collected in Table I Villi It had been 
argued in Ref. 86] that in the case of the noble met- 
als deviations from electron dynamics in a free gas of 
sp electrons due to the participation of d electrons in 
the screening of e-e interactions are of key importance in 
the determination of the inelastic lifetime of bulk elec- 
tronic states. Hence, Kliewer et al. [29J added this ef- 
fect to the calculated r e _ e foll owing the approach origi- 
nally suggested by Quinn 114j : they concluded that the 
screening of d electrons reduces the e-e scattering consid- 
erably, thus improving the agreement with experiment. 
Nevertheless, it was shown in Ref. (5(| that in the case 
of Shockley states, whose decay is dominated by intra- 
band transitions that are associated with very small val- 
ues of the momentum transfer, the screening of d elec- 
trons is expected to reduce the lifetime broadening only 
very slightly. Indeed, adding to the estimated Cu(100) 
Shockley e—e linewidth at T reported recently in Ref. [53| 
(with no (i-screening reduction) the e-ph linewidth of 7 
meV reported in Ref. [34J, the calculated total linewidth 
is found to be T ca i c = 26 meV in close agreement with the 
exprimcntally measured linewidth of 24 meV, as shown 
in Table IVIIII An extension of the approach reported in 
Ref. [53| to the case of the other noble metals Ag and 
Au should yield calculated linewidths that are closer to 
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experiment than those reported in Refs. [2!| and |l08| . 

The lifetime broadening of excited Shockley electrons 
beyond the f point (with k ^ and energies above the 
Fermi level - see Fig. [lj was studied with the ST M by 
Biirgi et al. [H on Cu(lll) and by Vitali et al. [Tl5| 
and Kliewcr et al. [116] on Ag(lll). The corresponding 
G°W° calculations that follow the I D sch eme of Chulkov 
et al. [4{| were reported in Refs. |ll7| and [ll5| | for 
Cu(lll) and Ag(lll), respectively, but now accounting 
for the potential variation parallel to the surface by intro- 
ducing not only a realistic effective mass for all bulk and 
surface states but also surface-state orbitals that change 
with k along the surface-state dispersion curve. More 
accurate calculations were later reported in the case of 
Cu(lll) [53| with the use of the k-dependent ID model 
potential of Eq. (|28p , showing that the inelastic lifetimes 
of excited Shockley electrons happen to be very sensitive 
to the actual shape of the surface-state single-particle or- 
bitals beyond the T point. A comparison between these 
more refined calculations and experiment demonstrated 
that there is close agreement at the surface-state band 



edge, i.e., at f , as shown in TableslVlland lVIlTl and there 
is also reasonable agreement at low excitation energies 
above the Fermi level. At energies where the surface-state 
band merges into the continuum of bulk states, however, 
the calculated lincwidths are found to be too low, which 
should be a signature of the need of a fully 3D description 
of the surface band structure. 
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